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The resonant-state expansion (RSE) Born approximation, a rigorous perturbative method de¬ 
veloped for electrodynamic and quantum mechanical open systems, is further developed to treat 
waveguides with a Sellmeier dispersion. For media that can be described by these types of disper¬ 
sion over the relevant frequency range, such as optical glass, I show that the perturbed RSE problem 
can be solved by diagonalizing a second-order eigenvalue problem. In the case of a single resonance 
at zero frequency, this is simplified to a generalized eigenvalue problem. Results are presented using 
analytically solvable planar waveguides and parameters of borosilicate BK7 glass, for a perturbation 
in the waveguide width. The efficiency of using either an exact dispersion over all frequencies or an 
approximate dispersion over a narrow frequency range is compared. I included a derivation of the 
RSE Born approximation for waveguides to make use of the resonances calculated by the RSE, an 
RSE extension of the well-known Born approximation. 

PACS numbers: 03.50.De, 42.25.-p, 03.65.Nk 


I. INTRODUCTION 

Fundamental to scattering theory, the Born approxi¬ 
mation consists of taking the incident field in place of 
the total field as the driving field at each point inside 
the scattering potential, it was first discovered by Born 
and presented in Ref. |l| . The Born approximation gave 
an expression for the differential scattering cross section 
in terms of the Fourier transform of the scattering po¬ 
tential. An important feature of this appearance of the 
Fourier transform is the availability of the inverse Fourier 
transform operation for the inverse scattering problem. 
The Born approximation is only valid for weak scatter- 
ers. In this paper I apply the RSE Born approxima¬ 
tion Ref. 00 to optical fibers or general open waveguide 
systems, which allows an arbitrary number of resonant- 
states (RSs) to be taken into account for scattering and 
transmission perturbative calculations. 

Optical fibers provide a well-controlled optical path 
which can carry light over potentially very long distances 
of several hundred to thousands of kilometers since the 
light is contained within the fiber by total internal re¬ 
flection. Fibers are clearly important for telecommuni¬ 
cation since the light may be modulated to carry infor¬ 
mation. It has recently been reported that a hollow-core 
photonic-band gap fiber yielded a record combination of 
low loss and wide bandwidth 0. Beyond telecommuni¬ 
cations waveguides are used in integrated optical circuits 
0, and terabit chip-to-chip interconnects [6]. 

Critical to understanding the response of a waveguide 
to being optically driven are its resonant states (RSs). 
The concept of RSs was first conceived and used by 
Gamow in 1928 in order to describe mathematically the 
process of radioactive decay, specifically the escape from 
the nuclear potential of an a particle by tunneling. Math¬ 
ematically this corresponded to solving Schrodinger’s 
equation for outgoing boundary conditions (BCs). These 
states have complex frequency w with negative imaginary 
part meaning their time dependence exp(— iuit) decays 


exponentially, thus giving an explanation for the expo¬ 
nential decay law of nuclear physics. The consequence of 
this exponential decay with time is that the further from 
the decaying system at a given instant of time, the greater 
the wave amplitude. An intuitive way of understanding 
this divergence of wave amplitude with distance is to no¬ 
tice that waves that are further away have left the system 
at an earlier time when less of the particle probability 
density had leaked out. There already exists numerical 
techniques for finding eigenmodes, such as the finite ele¬ 
ment method (FEM) and finite difference in time domain 
(FDTD) method to calculate resonances in open cavities. 
However, determining the effect of perturbations, which 
break the symmetry, presents a significant challenge as 
these popular computational techniques need large com¬ 
putational resources to model high-quality modes. Also 
these methods generate spurious solutions, which would 
damage the accuracy of the RSE Born approximation if 
included in the perturbation expansion basis. 

In order to calculate the resonances of open systems 
we have recently developed the resonant-state expansion. 
Such an approach was not previously available due to 
the lack of a normalization for resonant states. I derived 
the normalization of resonant states as a contribution 
to Ref. 0, and the generalization of the normalisation 
to dispersive media I derived straight forwardly at the 
time of my rigorous derivation. The problem of normal¬ 
ising resonant states stems from the fact that RSs with 
complex frequencies have wave functions which are expo¬ 
nentially growing in space away from system, making a 
volume integral over all space as would be done for a her- 
mitian system a meaningless exercise (see Appendix F). 
However the S -wave normalization was previously avail¬ 
able and analytically correct but numerically unstable as 
I showed analytically in Ref [3j]. 

So far the RSE has been applied to non dispersive sys¬ 
tems of different dimensionality IHH- However, almost 
all realistic systems have a relevant frequency dispersion 
of the refractive index. I have recently found 0,0 that 
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an Ohm’s law dispersion, i.e. a term in the susceptibility 
scaling at the inverse frequency, can be introduced to the 
RSE while keeping its linearity. We found in Ref. eh this 
dispersion can be a reasonable approximation for some 
materials over a limited range of wavelengths such as 
SHOTT BK7 glass over the optical range. In this work 
I generalize the RSE approach for waveguides detailed 
in Ref. [ 3 ] to systems constructed from dispersive media. 
Specifically I treat dispersion fully described by the Sell- 
meier equation, using a method similar to Muljarov et 
al in Ref. [13] for nano particles obeying Drude-Lorentz 
dispersion, however in this case generalised to inclined 
geometry. I also treat dispersion linear in wavelength 
squared, a generalisation of my work in Ref. id to in¬ 
clined geometry. I find the generalization of my Ohm’s 
law approach to inclined geometry to be greatly supe¬ 
rior to the generalization of the full dispersion treatment. 
This is most likely due to the unphysicality of the full dis¬ 
persion at high and low frequencies, which are beyond the 
fitting range of the dispersion models. In order to make 
use of the RSs, I derive the RSE Born approximation for 
waveguides Si, a theory for finding the field outside 
of a waveguide being internally or externally driven (see 
Appendix G). 

The paper is organized as follows, Sec. |TT| outlines the 
general recursive solution of Maxwell’s equation using 
Green’s function. Sec. IIIII develops the solution in Sec. El 
into a perturbation theory for waveguides with linear dis¬ 
persion in wavelength squared. Sec. |W] develops the solu¬ 
tion in Sec. El into a perturbation theory for waveguides 
with Sellmeier dispersion. Sec. I VII demonstrates and com¬ 
pares the two approaches for including dispersion into the 
RSE. The perturbation considered corresponds to a nar¬ 
rowing of the waveguide and is discussed in Sec. IVI Cl 
In Sec. IVI Dl I show the results using the simple disper¬ 
sive approximation BK7. In Sec. IVI El I use the full Sell¬ 
meier dispersion of BK7. Finally, we discuss the compar¬ 
ison and performance of the two methods in Sec. IVI Fl in 
Sec. Ivm I derive the equations for the RSE Born approx¬ 
imation for planar waveguides. In Appendix C, D, E, and 
F I generalize my part of the proof of the normalization 
of RSs to open waveguide systems. 


II. RESONANT STATE EXPANSION FOR 
WAVEGUIDES AND NON-NORMAL 
INCIDENCE 


In this section I develop the general recursive solution 
to the problem of calculating the resonant states of a per¬ 
turbed waveguide. The recursive solution requires the 
Green’s function (GF) of the unperturbed waveguide, 
and a perturbation which is retaining the translational 
invariance along the waveguide. 

I first consider a waveguide of thickness 2 a in vacuum 
with translational invariance in one direction, having the 


dielectric constant 

Mr) 


for |r| < a , 
1 for |r| > a , 


(1) 


where e u is the frequency dependent relative permittivity 
(RP) of the wave guide with the angular frequency w. 
r are the coordinates normal to the waveguide, which 
for a cylindrical waveguide is given in polar coordinates 
r = (p , 9) and for a planar waveguide r is Cartesian x. I 
assume a relative permeability of /i = 1 throughout this 
work. The electric field E satisfies Maxwell’s equation, 


E(r,f) = 0, 


( 2 ) 


For cylindrical and planar waveguides, due to transla¬ 
tional invariance in the z direction I first assume then 
prove that E can be factorised as follows 

E(r,i) = e^-^E(r), (3) 


in which p is the wave vector along the translationally 
invariant direction. For the component E(r) of the elec¬ 
tric field, Eq. (O transforms to a one-dimensional (ID) 
wave equation in the planar case or a 2D wave-equation 
for the cylindrical waveguide to 


-L + e w (r) 


LU 


E(r) =0. 


( 4 ) 


Here L is a linear operator not dependent on z or e u (r). 
The form of L is derived in Appendix C, the derivation 
shows that the trial solution Eq. ([3]) is correct since it 
demonstrates that E can be factorised by separation of 
variables. In principle L could be any linear operator 
independent of both z and e w (r), and hence the treat¬ 
ment of waveguides in this paper is readily applicable to 
quantum mechanics or acoustics. 

The non-normal incidence, characterized by p ^ 0, is 
treated here. The previously used spectral representa¬ 
tion of the GF in the frequency domain contains a cut 
for p ^ 0, which can be removed by mapping the problem 
onto the complex normal wave-vector space fc, as demon¬ 
strated in Ref. 0 - The relation between k, p, and ui is 
defined by us to be u 2 /c 2 = k 2 +p 2 , k and p are or¬ 
thogonal by Pythagorean identity. I follow also here the 
approach used in Ref. 3 and formulate the RSE in the 
complex fc-plane, for which the spectral representation 
of the GF of an infinite planar system with an in-plane 
momentum p / 0 written in the spectral form 


Gi 1) (r,f') = E 


E n (r) ® E n (r') 

2 k n (k - k n ) 


( 5 ) 


where E n (r) is the electric field of a RS, defined as an 
eigensolution of Eq. (|4|) with an arbitrary profile of i u (r) 
within a region |r| < a, satisfying the outgoing wave 
boundary conditions at infinity. The dyadic product of 
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the E„(r) fields is denoted by <g). The in-plane eigen- 
frequency k n are the poles in Gfc(r, r'). By definition 
w 2 /c 2 = + p 2 where w 2 is the eigenfrequency corre¬ 

sponding to the eigenstate E n (r). For a study of the 
derivation of such spectral representations see Appendix 
D. 

The GF satisfies the equation 

[-L + £w(r )(k 2 + p 2 )] Gfc(r, r') = i<5(r - r') (6) 

In the present work, we consider the refractive profile 
(RP) having the property 

lim e w (r) = e(r), (7) 

uj —>-oo 


where e(f) if the frequency independent part of e w (r). 
If we assume i u (r) to be discontinuous through the 
boundary of the system then by substituting Eq. (0 into 
Eq. 0. convoluting with a finite function, and letting 
h —> oo, we obtain the sum rule (see Appendix E) 


E 


E„(r) <g> En(r') 


k 


n 


= o, 


( 8 ) 


which allows us to re-write the Green’s function a second 
way as 

p( 2) / f f /\ V- E n (r) ® E n (r') f 1 _ k 

k 1 ’ > 07 . k — k k 2 + n 2 ’ 

(9) 

I now consider an arbitrary perturbation Ae u (r) of the 
dielectric constant inside the layer |r| < a. I use the 
Eq. (0 to solve the perturbed problem, 

£ v (y) = J G k (r,r')AUr')SAr')dr', ( 10 ) 

which is the solution to the equation 

L£ v (r) = [e u (r) + Ae u (r)] £„(r) (11) 

c A 

Note that the perturbed modes £„(r) satisfy Eq. (J4|) 
with e w (r) replaced by £ u (r) + A£^(r) and the outgoing 
boundary conditions with ui n replaced by u>. I show the 
solution of this equation perturbatively for two different 
types of dispersion in the following two sections. 


III. SINGLE SELLMEIER RESONANCE AT 
ZERO FREQUENCY 


of the Green’s function, similar to Ref. [L3J , Eq. 0 and 
Eq. 0 for the different components of the perturbation 
in Eq. (flOl) , 


£{v) 


— J G ( k \r,r')Aa-{r')£{r')dr' 

2 J Gfc 2) (f,r')A£(r ■')£{r')dr\ (13) 


which results in the following relationship between un¬ 
perturbed and perturbed modes 


En(f) 


S(r) = -(k 2 +p 2 )J2 2k k _ k k2 o 

^jrV'Yl |_ E •' J Tl »V | JJ 

x f E n (r , )Ae(r , )£(r' ')dr 
J —a 


E„(r) 


E 2 k n (k - k n ) 

x f E n (r')A(x(r')£(r')cZr. (14) 

J—a 


The perturbed mode £(r) satisfies Eq. 0 with £ u (r) re¬ 
placed by £ u (r) + Ae u (r) and outgoing boundary con¬ 
ditions with the wave vector k. In the interior region 
|r| < a which contains the perturbation, the perturbed 
RSs £ v of wavenumbers x v can be expanded into the un¬ 
perturbed ones, exploiting the completeness of the latter 
which follows from Eq. (IA6I) (see also Appendix E), 


£/(r) = E 6 ^ En (r) • (15) 

n 


Substituting this expansion into Eq. (fl3l) and equating 
coefficients at the same basis functions E n results in the 
matrix equation 






( P 4 ^ifk n )V n m 4 ” 
k n ) 


where 


and 


Vnm = J E n (r) • A£(r)E m (r) dr 


S nm = / E n (r) • Aff(f)E m (r) dr. 


(16) 


(17) 


(18) 


With the substitution c nv = b mk /k^, Eq. (fTTH) can be 
rewritten as 


For a Sellmeier dispersion with a single resonance at 
zero frequency (SRZ), I write the perturbation in Eq. (flOl) 
as 

c 2 Ao’(r) 

AM*) = A£(r) + C (12) 

UJ A 



again present only inside the layer |r| < a. It is then care should be taken to take the sign of \fk, „ consistently 
possible to linearise the RSE by using the different forms between matrix elements. Eq. (1191) is linear in n v and 
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can be solved by libraries for generalized linear matrix 
eigenvalue problems. In the absence of dispersion, S nrn = 
0, and Eq. (1191) reverts back to the expression for non- 
dispersive waveguides [7]. In the absence of p, pa = 0, 
we see that Eq. m reverts to an expression for dispersive 
perturbation to nano-particles a 


IV. SELLMEIER DISPERSION 


In this section I develop the existing RSE for wave 
guides into a perturbation theory for waveguides with 
Sellmcier dispersion, following a similar approach to 
Ref. [Till] . The Sellmeier dispersion is the sum of a set 
of Lorentzians in frequency with poles qj on the real 
axis. The starting point for this derivation is Eq. m, 
the recursive solution for the perturbed problem derived 

in Sec. HIl 

In this section, we consider the RP 


Mr) + Y1 

3 


c 2 M r) 


( 20 ) 


and perturbation, 


Ae w (r) 


Ae(r) + Y^ 

3 


c 2 A&j(r) 

u> 2 - nj 


( 21 ) 


with j numbering the resonances at frequencies Qj hav¬ 
ing oscillator strengths &j. We introduce an effective 
resonance wave vector qj as q 2 = ft 2 /c 2 — p 2 and re-write 
Eq. m as, 


AMr) = A£(r) + ^ ^ J ^ 


( 22 ) 


= Ae(f)+^] 


( Affj(f) + Affj(r) 


\2k(k — qj) 2 k(k + qj) J 


In the Appendix A I use my method of [2], which was 
first adapted for full dispersion of 3D nano-particles by 
E. A. Muljarov in Ref. 13|, to derive for non-vanishing 
<r j the sum rule 


E n (r) ® E n (r / ) 

n 2/.' n ( k n dz qj) 


(23) 


which allows us to write the Green’s function as 


p(±j),. .A _ E n(r) ® E n (r')(fc ± qj ) 
k ’ M 2 k n (k — k n )(k n ± qj) 


which has the useful (fc ± qj) in the numerator for 
reducing the order of the perturbation matrix prob¬ 
lem through cancellation with perturbation denomina¬ 
tors which would otherwise have lead to high order poly¬ 
nomial eigen problems. 


We now consider a perturbation Ae u (t) of the RP in¬ 
side the layer |r| < a. Similar to the previous section we 
use the Eq. © to solve the perturbed problem, 


£(£■) =-- / Gfc(r, r')Ae ul (r')£(r')dr ', 


(25) 


where we take Ae w (r) of the same form as Eq. m- Us¬ 
ing both forms of the Green’s function HI, Eq. © and 
Eq. (1)H1) in Eq. (05l) . yields 


£(r) = - 


gJ, 2) (r, r')Ae(r')£ (r')dr' 


(26) 


UJ V "S 

±,j • 




Ad-^P) 

2 k(k ± qj) 


£(r')dr'. 


which results in the following relationship between un¬ 
perturbed and perturbed modes 


£{r) 


-(* 2 +«•’)£ 

n 



l 

k - kn 


k 

k 2 + p 2 


x J E, n (r')Ae(r')£(r')dr' 

- (k 2 + v 2 ) V E n (r)(fc ± qj) 
k ~^j 2k n (k n ±qj){k - k n ) 


x 


E„(r') 


Ad-j(r') 

2k(k ± qj) 


£(r')dr l . 


(27) 


The perturbed mode £(r) satisfies Eq. Q with e w (r) re¬ 
placed by £cu(r) + Ae w (r) and outgoing boundary con¬ 
ditions with the wave vector k. In the interior region 
|r| < a which contains the perturbation, the perturbed 
RSs £ v of wavenumbers x v can be expanded into the un¬ 
perturbed ones, exploiting the completeness of the latter 
which follows from Eq. (IA6I) (see also Appendix E), 


£„(r) = ^2b nu E n (r). 


(28) 


Substituting this expansion into Eq. (l27l) and equating 
coefficients at the same basis functions E n results in the 
matrix equation 


brtis 


where 


and 


p 2 + x v k n 
2 k n (^K v k n ) 

(xl+p 2 ) 

2>f v (yf v - k n ) ^ k 2 - q 

v ' m,j J 


^ ^ Vnmbn 

m 

_ A {j) 

V ^ -g* -nm , 

—- b 


2 u mv i 


V nm = J E n (r) • A£(r)E m (r) dr 


(29) 


(30) 


A nm = / E „(r) • Affj(r)E m (r) dr (31) 
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is the matrix of the perturbation in the basis of unper- RSE : {k n } —> (40) 

turbed RSs. Introducing the abbreviation 


t O') 


Unrr). 


j 


~ Qj) 


(32) 


I arrive at 





nm 




= o, 


(33) 


which is of second order in x u . We can write this ma¬ 
trix problem compactly as Q(>r„)b„ = 0 with Q(>r„) = 
>f 2 M+xyC+K. To solve_this matrix problem for a basis 
of size N, I follow Ref. 0 and use the first companion 
linearization, defining 


R(x„) = x, 


M 0 

0 I 


C K 

-I 0 


(34) 


where I is the N-by-N identity matrix, with the corre¬ 
sponding vector 

z = (35) 

We solve R(x„)z = 0 using the generalized eigenvalue 
solver from the numerical algorithms group (NAG) C-|—b 
library. We then take the first N components of z as the 
eigenvector b„. 


V. POSSIBLE EXPLANATION FOR THE 
DIFFERENCE IN ORDER BETWEEN THE TWO 
METHODS 


RSE : {E n } -> {£„} (41) 

where [e^o; 2 ]^ poles ^ and [Ae{,w 2 ]^ poles ^ have poles in the 
same position in frequency space while [e w a; 2 ]^ res ^ and 
[e^w 2 + Ae w w 2 ]( reg ) are regular functions. 

The mapping in Eq. (138)) is between two spaces which 
are topologically equivalent to non-dispersive spaces. 
Therefore, the mapping problem is mathematically 
equivalent to non-dispersive RSE perturbation theory 
and so the order of the eigenvalue problem is not in¬ 
creased upon that of the non-dispersive problem. 

However, in the case of mapping by Eq. (l39l) 
[e^w 2 ]( poles ) contains poles in frequency space and so 
cannot be continuously deformed into [E u] tx 2 ]( iee ' > hence 
the map given by Eq. (1551) is not equivalent to a 
non-dispersive mapping and so the order of the RSE 
eigenvalue problem is increased upon that of the non- 
dispersive problem. 

Formally A and B are topologically equivalent if there 
is a homeomorphism mapping orbits of A to orbits of 
B homeomorphically, and preserving orientation of the 
orbits. 


VI. APPLICATION TO A PLANAR 
WAVEGUIDE 

In this section I discuss the application of the Sellmeier 
waveguide RSE to an effectively ID planar waveguide 
system, translationally invariant in the Cartesian % and y 
directions, described by a scalar RP, i.e., i u (x) = 1e u (x), 
A i UJ (x) = lAe u (x). As unperturbed system I use a ho¬ 
mogeneous planar wave guide of half width a, so that 


Here I go beyond the algebraic mathematics presented 
so far and discuss the fundamental topological differences 
between the two RSE approaches that I have developed. 

The RSE perturbation theory presented here is an in¬ 
jective continuous mapping function except in the neigh¬ 
bourhood of the finite number of poles such that, 

e^ 2 = [e u uj 2 } ( reg ) + ^ [e£ w 2 ] (poles) (36) 

3 


( for \x\ < a, 
1 1 elsewhere. 


(42) 


A. Unperturbed resonant states 

The solutions of Eq. JI]), which satisfy the outgoing- 
wave boundary conditions in TE polarization take the 
form Q 


is mapped to 

£ u w 2 + Ae^w 2 = [e^uj 2 + Ae u w 2 ] (reg) (37) 

+ ^[e> 2 + A££a, 2 ]( poles \ 

3 

in the following way 

RSE : [e w w 2 ] (reg) -► [ e + A£^ 2 ] (res) (38) 


f (-1 ) n A n e- iknX , x < -a, 

E n {x) = l B n [e iqnX + (—l)"e - * 9nX ], \x\ < a, (43) 

[ A n e lknX , x ^ a, 

where the eigenvalues k n satisfy the secular equation 

(k n - q n ) e iqna + (-1)” (k n + q n ) e~ iq ~ a = 0 , (44) 

with 


RSE : [e> 2 ]( poles ) [e^w 2 + Ae> 2 ]( poles) (39) 


q n = V e^k 2 + (e„ - 1 )p 2 . 


(45) 
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I use here an integer index n which takes even (odd) 
values for symmetric (anti-symmetric) RSs, respectively. 
The normalization constants A n and B n are found from 
the continuity of E n across the boundaries and the nor¬ 
malization condition found in Appendix F. 

In order to arrive at the normalization condition I con¬ 
sider an effectively 2D system and I take E(r, k) as being 
the analytic continuation into k space of E n (r) and A 
being an arbitrary cross-section of the translational in¬ 
variant waveguide to arrive at (see Appendix F) 


1 + Sk n ,o = / E n (r) 


d(w 2 e^(r)) 


E n (r)dr 


+ lim 

k^k n 


<9(u/ 2 ) 

E n (r)LE(r,k)-E(r,k)TE n (r 


h.2 _ U2 

lx 


(46) 

dr 


outside the system in free space Maxwell’s equation sim¬ 
plifies (see Appendixes F and G), therefore by extending 
A outside the system and denoting its circumference La 
I can write, 


l + <5fc„,o — / E n (r) 


d(uj 2 £ u ( r)) 


<9(w 2 ) 


E n (r)dr(47) 


+ lim t E "~ VE - E . VE “dL. 


k—tk n 


U2 _ U2 

lx A/„ 


I made the necessary assumption that is a (real) sym¬ 
metric matrix or a scalar so that (defined to be in this 
case) E ■ e w „E n = E ra • e^„E and non dispersive at high 
frequencies (see Appendix E and F). 

Various schemes exist to evaluate the line integral limit 
in Eq. ( 1571 ) such as analytic methods in Ref.[1, [ll[ or 
numerically extending the surface into a non-reflecting, 
absorbing, perfectly matched layer where it vanishes. 
Hence, I derive from Eq. (1471) the relevant normalization 
condition for the planar waveguide systems, which I will 
use for the numerical demonstration, 


5(w 2 £^(a;)) 


- a d{^) 

E n (-a)E n (-a) + E n (a)E n (a) 
i2k n 


E n (x)E n (x) dx 

= l + 4 n ,o.(48) 


Here in Eq. (l48l) the first integral is taken over an arbi¬ 
trary simply connected line enclosing the inhomogeneity 
of the system and the center of the coordinates used, and 
the second term is evaluated at its end points. 

The coefficients in Eq. ( 1551 ) take the form 


H) n 

\J 4(ae w + ip 2 /(k n u} 2 /c 2 )) + M n r] ’ 


(49) 


with 


An — Bn 


(—l) n e i9 " a 


1 <9(w 2 e<J 
11 ~ ^ d(w 2 ) 


- 1 , 


(50) 


(51) 


and 


M n = 


(-l) n 2sin(2 q n a) 
-b 4a 

Qn 


(52) 


where w 2 /c 2 = k 2 + p 2 ■ 

Around each pole of e u the secular equation Eq. (l44l) 
has a countable infinite number of solutions, each one 
creating a RS. An analytic approximation to these solu¬ 
tions is given in Appendix B, which are used as starting 
values for the numerical solution of Eq. (l44l) to determine 
these RSs closely spaced in the complex frequency plane. 


B. BK7 Sellmeier dispersion 


In our numerical examples I use dispersion parameters 
describing the common borosilicate glass SCHOTT BK7. 
Its refractive index n r in the optical frequency range is 
well described by the Sellmeier expression 


€(jJ 


1 + 


1.03961212A 2 
A 2 - 6000.69867 nm 2 
0.231792344A 2 
A 2 - 20017.9144 nm 2 + 


(53) 

1.01046945A 2 
A 2 - 103.560653 ^m 2 ' 


6839.577136 

(aw/c) 2 - 6578.970180 ( ' 

457.1302872 0.3852016551 

_ ( auj/c) 2 - 1972.154382 _ (aw/c) 2 - 0.3812105899 

The resulting RP is real and is shown in Fig.|T] This dis¬ 
persion is of the form Eq. (l20l) and can be treated using 
the quadratic matrix equation Eq. ( 1551 ) . In order to com¬ 
pare this result with the one of the linear matrix equation 
Eq. (fl9l) . I have fitted by e u over the wavelength range 
from 1.25 //m to 1.75 pm with the form of Eq. m, yield¬ 
ing 

£u> = 2.28239 - 0.01262 X 2 /pm 2 (55) 

The deviations of this fit over the fitted range are around 
10 -4 , as shown in Fig.|T] 


C. Unperturbed and perturbed system 

The unperturbed system I consider in this work as ex¬ 
ample is a planar waveguide of width 2a = 2 pm. I con¬ 
sider a size perturbation, narrowing the waveguide by 10 
percent, as shown in Fig. [2] For the in-plane wave vector 
component I have chosen pa = 5. I choose as basis of size 
N for the RSE all RSs with 

IfcnV^wTl < fcmax(A0 , (56) 

i.e., with a wave vector in the medium below a suitably 
chosen maximum k mSbX {N), this follows the approach of 
Ref. HI. 
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FIG. 1: RP e u of SCHOTT BK7 glass given by the Sellmeier 
expression Eq. (1531) (solid line), and fitted using a single reso¬ 
nance at zero frequency (dashed line). The absolute residual 
of the fit is given as dotted line. 



wavelength X (pm) 



<| Q-' _I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_ 

3.0 3.5 4.0 4.5 5.0 5.5 6.0 

RS frequency Re(aro/c) 


FIG. 2: Cross-section of the perturbed and unperturbed wave 
guide used in Sec. IVI Dl and Sec. IVI El shown at a particular 
frequency. 


The motivation for the choice of basis selection criteria 
is the analogy between the RSE and a Fourier expansion. 
As I increase k n y/e UJn , I increase the number of oscilla¬ 
tions in the field as can be seen from Eq. (l43l) . These in¬ 
creasingly oscillating fields go into the basis and similar 
to Fourier series, increase the resolution of the composite 
field generated by the expansion. 


D. Single Sellmeier resonance at zero frequency 

Here I show results of the SRZ yielding the linear eigen¬ 
value problem Eq. m- For the unperturbed and per¬ 
turbed system I use a RP given by the SRZ Eq. (l55l) . 
The refractive index of the unperturbed and perturbed 
system is compared in Fig. [3] with the Sellmeier expres¬ 
sion Eq. (l53l) . I can see that for the chosen waveg¬ 
uide width, the fitted wavelength range corresponds to 
3.6 < aui/c < 5. The RS frequencies are given in Fig. [3Kb) 
for the unperturbed system uj n = cy^ + p 2 and for the 
perturbed system vj v = C\J yfy + p 2 . To compare 


FIG. 3: SRZ RSE results for a thickness perturbation of 
a planar waveguide as function of the real part of ui. (a) 
Refractive index of the unperturbed and perturbed medium 
(dashed line) described by the SRZ Eq. (I55|) and of the BK7 
Sellmeier dispersion (solid line) Eq. (1531) . (b) RS frequencies 
for N = 200. Shown are exact unperturbed (open circles), ex¬ 
act perturbed (open squares), and RSE perturbed (crosses) 
data, (c) Relative error A of the RSE perturbed RS frequen¬ 
cies for N = 50,100, 200. Closed symbols give relative error to 
the SRZ Eq. (1551) RS eigen-frequencies as shown. Open sym¬ 
bols give relative error to the Sellmeier dispersion RS eigen- 
frequencies as shown. I show in (b) states with both incom¬ 
ing and outgoing boundary conditions, i.e. having positive or 
negative imaginary parts of complex resonant frequency, only 
states with outgoing boundary conditions go into the basis. 

calculated using the RSE with the exact result 
obtained by the secular Eq. (l44l) , I define the relative er¬ 
ror A = „/vo l 'v Xaci ' 1 — 11, and give the resulting values 

in Fig.[3Kc). I find that as I increase N , A decreases 
proportional to N ~ 3 , similar to the findings for the non- 
dispersive RSE P~[lOl|. and values in the 10 -6 range are 
reached for N = 200. This is actually smaller than the 
relative error due to the SRZ approximation of the Sell¬ 
meier dispersion (see Fig. [I]). 

The evolution of the perturbed RSs wavenumbers w v 



























in-plane wavevector pa in-plane wavevector pa 


FIG. 4: SRZ RSE results for a thickness perturbation of a 
planar waveguide and N = 50. Shown are K(u)„) as function 
of the in-plane wave vector p with given by the color 

according to the scale shown. 


FIG. 5: SRZ RSE results for a thickness perturbation of a 
planar waveguide and N = 50. Shown are as function 

of the in-plane wave vector p with the relative error A given 
by the color according to the scale shown. 


with the in-plane wave vector p is shown in Fig. [3] I can 
distinguish j7J the waveguide (WG) and anti-waveguide 
(AWG) modes, which have real m v , and the Fabry-Perot 
(FP) modes which have a finite imaginary part represent¬ 
ing their losses. Increasing p, the FP modes split into 
WG and AWG mode at the bifurcation point at which 
pc = vj v , at which k = 0 i.e. at grazing incidence of the 
external field. 

The relative error A of w v is given in Fig.[5]as function 
of p. I can see that A has generally a weak dependence 
on p, except close to the bifurcation points, where the er¬ 
ror of the FP and AWG mode is significantly increased. 
Under closer examination I see further smaller disconti¬ 
nuities as function of p, which correspond to a change of 
the basis states included according to Eq. (l56l) for fixed 
N = 50. These results indicates that the RSE is able to 
reproduce all relevant RSs with a good accuracy. 


E. Sellmeier Dispersion 

Here I show results of the Sellmeier RSE yielding the 
quadratic eigenvalue problem Eq. ( 1351 ) . 

The RS frequencies are given in Fig.^b) for the un¬ 
perturbed system uj n = C\Jk ^ + p 2 and for the perturbed 
system = c*J + p 2 . The relative error A is given 
in Fig. He). Also here I find that as we increase N, A 
decreases proportional to TV -3 , and values in the 10~' 5 
range are reached for N = 800. 

The evolution of the perturbed RSs wavenumbers w v 
with the in-plane wave vector p is shown in Fig. [3 in¬ 
cluding the relative error A, and in Fig. [I] including the 
imaginary part 3(t*7„). 

I can see that A has generally a weak dependence on p, 


except for the bifurcation points, at which a Fabry-Perot 
mode splits into a waveguide and anti-waveguide mode. 

These results indicate that with full Sellmeier disper¬ 
sion the RSE is able to reproduce all relevant RSs with 
a good accuracy. 

F. Performance Comparison 

Here I compare the computational complexity of the 
SRZ RSE and the Sellmeier RSE. The corresponding 
eigenvalue problems Eq. (USD and Eq. (1331) show that for 
the same number of basis states N, the SRZ RSE is a gen¬ 
eralized eigenvalue problem of size N x N, while the Sell¬ 
meier RSE has a size 2 N x 2 N due to the quadratic nature 
of Eq. (|33l) . Also I see from Fig.[8]that for pa = 5 the reso¬ 
nances closely associated to the Sellmeier qj poles, those 
eigen-frequencies found using the solutions to Eq. (IB4I) 
as a starting point for the Newton-Raplison search, con¬ 
tribute approximately 10% to the basis size thereby in¬ 
creasing the numerical complexity and reducing the effi¬ 
ciency. 

In Fig. [9] I show a comparison between the average 
relative error in the perturbed resonances for pa = 5 
found in the range 3 < aui < 5 versus the number of 
seconds required to diagonalise the perturbation eigen¬ 
value problems using an Intel Core 2 Duo, 6M Cache, 
3.16GHz, 1333MHz FSB processor connected to 4GB of 
RAM and using the NAG generalized eigenvalue problem 
solver software. 

For the case of SRZ dispersion the relative error reaches 
the limit of about 5 x 1CU 5 due to the RP approximation 
at A « 70. Once the perturbation method using the full 
Sellmeier dispersion reaches this relative error I see that 
upon adding further basis states more noise in the relative 
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3456789 10 

RS frequency Re(ato/c) 


FIG. 6: Sellmeier RSE results for a thickness perturbation 
of a planar waveguide, (a) refractive index of the Sellmeier 
dispersion given by Eq. (15311 . (b) RS frequencies for N = 800 
basis states. Shown are exact unperturbed (open circles), ex¬ 
act perturbed (open squares), and RSE perturbed (crosses) 
data, (c) relative error A of the RSE perturbed RS frequen¬ 
cies V3 U for N = 50,100, 200, 400, 800 basis states, as labeled. 
I show in (b) states with both incoming and outgoing bound¬ 
ary conditions, i.e. having positive or negative imaginary 
parts of complex resonant frequency, only states with out¬ 
going boundary conditions go into the basis. 


error plot is generated and that the average relative error 
changes little. Fig.[9]suggests that the SRZ RSE is several 
orders of magnitude more efficient than the full Sellmeier 
RSE when considering the resonances over a small range 
of frequencies such as the small range used for optical 
communications. 



in-plane wavevector pa 


FIG. 7: Sellmeier RSE results for a thickness perturbation 
of a planar waveguide and N = 200. Shown are 5 H(zj 7,/) as 
function of the in-plane wave vector p with the relative error 
A given by the colour according to the scale shown. 



FIG. 8: Graph showing quantitatively how the number of 
poles in the basis associated closely with the poles in the Sell¬ 
meier equation for pa = 5, those eigenfrequencies found using 
the solutions to Eq. (IB41) as a starting point for the Newton- 
Raphson search, increase with basis size N when the basis 
is selected according to Eq. (15611 . For reference I also show 
the number of ordinary poles, those found with the Newton- 
Raphson method without requiring the solutions of Eq. (lB4l) 
as a starting point for their search and discovery, in the basis 
versus basis size N. 


VII. RSE BORN APPROXIMATION FOR 
PLANAR WAVEGUIDES 


In this section I derive the RSE Born approximation 
for planar waveguide, the equivalent theory for effectively 
2D systems is given in Appendix G. 

In effectively one dimension Eq. @ becomes in free 


space (taking c = 1), 



1 d 2 ' 
c 2 dt 2 


e p E{ x ) = o, 


However, V 2 0 p = — p 2 9 p , therefore 


dx 2 


k 2 


E{x) = 0, 


(57) 


(58) 
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FIG. 9: Comparison of the average relative error of the wave 
guide modes occurring in the range 3 < au) < 5 calculated 
using the two different perturbation schemes in this paper 
versus time required to diagonalise the perturbation matrix 
using an Intel Core 2 Duo, 6M Cache, 3.16GHz, 1333MHz 
FSB processor connected to 4GB of RAM and using the NAG 
generalised eigenvalue problem solver software 


Hence the free space GF equation is, 


- 

dx 2 


G[ s (x, x') = S(x — x') 


(59) 


which has the solution, 


G{ s {x,x') = -- 


0 ik\x—x'\ 


2 ik 


(60) 


The systems associated with G{ s and Gk of Eq. (JG)) are re¬ 
lated by the Dyson Equations perturbingback and forth 
with Ae ul (x) = e UJ (x) — 1 similar to Ref. [4], 


G k {x,x") 


G{ s (x,x") (61) 

—w 2 J G{ s (x, x"')A,£ u (x'")G k (x'", x”)dx" , 


Hence, in one dimension the RSE Born approximation 
can be greatly simplified by using Eq. (l60l) and the spec¬ 
tral GF Eq. (1E7I) in Eq. (1631) to arrive at, 

p ik(x" — x) (.2 p ik{x" — x) pa 

w 4 e ifc(x"-x) A n (x)A n (x") 

4 k 2 2s 2 k(k-kn) ' 1 ’ 

n v y 

where A n is defined as the Fourier transform, 

A n (x)= f e ik&x ’Ae u {x')E n (x')dx' (65) 

J —a 

Interestingly in one dimension I do not require the 
far field approximation to make the simplification of the 
Green’s function required to bring the RSE Born approx¬ 
imation to the form of Eq. (IG9I) . Hence, in one dimension 
the RSE Born Approximation is valid everywhere outside 
of the slab and not just in the far field. I note that fast 
Fourier transforms are available for use upon Eq. (1651) . 

I demonstrate the computational accuracy of the RSE 
Born approximation in Fig. [10] where I calculate the 
transmission defined as 

T(k,x') = \2kG(x',-ci-k)\ 2 (66) 

For comparison the analytic GF is found by solving 
Maxwell’s wave equation in one dimension with a source 
of plane waves while making use of Maxwell’s bound¬ 
ary conditions. The system treated is the unperturbed 
system in Sec. I VI Dl with pa = 5. From Fig.[TO] we can 
see that unlike the standard Born approximation the 
RSE Born approximation is valid over an arbitrarily wide 
range of k depending only on the basis size N used. Fur¬ 
thermore we see that as the basis size increases the RSE 
Born approximation converges to the exact solution. The 
absolute error in the RSE Born approximation is approx¬ 
imately reduced by an order of magnitude each time the 
basis size is doubled. Absolute errors of10 -8 — 10~ 5 are 
seen in the k range shown for basis size N = 401. 


G k (x"',x") = G{ s (x"',x") (62) 

-w 2 J Gk{x"', x')A£ U] (x')GI s (x ', x")dx' , 

Combining Eq. (1CT1) and Eq. (IS^l) it is obtained similar 
to Ref. [2[ 

G k (x,x") = G{ s (x, x") 

- W 2 J G f k s (x, x^Asuix^Gl 3 (x 1 , x")dx' 

+ w 4 J J G f k s (x,x')A£ ul (x')G k (x',x r ") 

x A£ uj (x"')gI s ( x" r ,x")dx"'dx' . (63) 


VIII. SUMMARY 

In this work I have extended the RSE to media having 
a simple dispersion linear in wavelength squared. This 
dispersion has a single pole at zero frequency and does 
not introduce an additional dynamic degree of freedom as 
it would be the case for a more general Sellmeier model 
of the material response. This property allows to keep 
the simplicity of the RSE formulation, therefore retain¬ 
ing the advantage of the RSE in computational efficiency 
discussed in Ref.H. 

Furthermore, in this work I have also extended the 
RSE to waveguides obeying the Sellmeier equation for 
glasses. In order to do this I have reduced the order 























11 



frequency ( aw/c) 


FIG. 10: I give transmission results for the unperturbed sys¬ 
tem in Sec. IVI Dl with pa = 5. I have used analytic modes 
as RSs for the RSE Born approximation, (a) Exact trans¬ 
mission as a function of frequency as defined by Eq. (16611 (b) 
Absolute error in transmission calculated using the analytic 
form of T(k,x') between x = —a and x = a versus fre¬ 
quency uj as comparison for the RSE Born approximation. 
Here N = 51,101, 201, 401 as labeled. 


Appendix A: Sum rule 

Here I use my approach of Ref.Q, which was first 
adapted for full dispersion of 3D nano-particles by E. A. 
Muljarov in Ref. [131] . to derive sum rules and complete¬ 
ness of modes for the waveguide GF. 

By substituting Eq. (0 into Eq. 0 and using Eq. 0 I 
have, 


- id(r-r') = 

n 


c 2 2k n (k - k n ) 


E n (r) <g> E n (r'). 


(Al) 


I now consider the pole in e u at k = 

= Qj-(r) 

“ 2 k(k±qj) 


(A2) 


where I have been able to ignored other terms in the 
dispersion since they are constant in the limit k —> =p qj. 
Using w 2 /c 2 = k 2 +p 2 and w 2 /c 2 = k 2 +p 2 and Eq. (IA2I) 
in Eq. (1 Al [1 . I find 


((U 2 ^ - CJ 2 £uj) 

c 2 2k„(k - k n ) 

_ /.x__ 

<Tj r 8 k 2 k(k - k n )(k n ± qj) 
2*j(j) ( P 2 kn k n q 0 k \ 

( k±qj ) \ 8 k 2 k{k-k n ) 8 k 2 k(k n ±qj)J 


Then convoluting Eq. ED with arbitrary finite function 
of the eigenvalue problem to second order by using sum an d taking the limit k -A =F Qj, I see that the second term 
rules as was done in [e| . is diverging unless 


To make use of the RSs generated by the RSE I derived 
the RSE Born approximation for effectively ID and 2D 
waveguides. 

My conclusion is that the simple dispersion treatment 
is more efficient than using the full Sellmeier dispersion. 
This is because the dispersion only has to be correctly 
reproduced over a narrow part of the optical frequency 
range and therefore it is inefficient to include the full and 
at some frequencies unphysical Sellmeier dispersion. 


E 


E n (r) ( 8 ) E n (r') 
k n (k n ± qj) 


= 0. 


since when k —> qj 


(A4) 


( p 2 k n 


± 


k n qjk 


\8 k 2 k(k — k n ) 8k 2 k(k n ±qj) 

= ± 


8 k n {k n ± qj)qj 


(AS) 
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Appendix B: Analytic solution of secular equation 
close to RP poles 

Here I calculate an analytic approximation for the so¬ 
lutions of the secular equation Eq. CH) close to the poles 
of the RP, following the approach of Ref. jlSj]. In the limit 
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—► oo I find from Eq. (1451) that q n -A yfcZiiOnl c, so that 
Eq. (TH1) is given by —e l9 " a + (—l) n e -l9n “ = 0, and thus 


sin (y/e^uj n a/c) = 0 for even n 
cos {y/e^uj n a/c) = 0 for odd n 


(Bl) 


In the limit k n —> ±qj to the poles of the RP given in 
Eq. (1221) I find 


u 2 n P n +p 2 

u c 2 J 2k n (k n =F Qj) 

and solutions of Eq. m are given by 

7 no 2 = (t rn/2) 2 


(B2) 


(B3) 


The solutions of Eq. (IB2I) with k n close to ±qj are given 
by the analytic solution of the quadratic formula 


(tjj - 2y n )kl ± 2y n qjk n + <Jjp 2 = 0 (B4) 

I use the solutions of Eq. (1B4I) as starting points for the 
Newton-Raphson search for RSs. 

The RSs of the unperturbed system (see Sec. I VI Cl) 
with Sellmeier dispersion Eq. (1531) are given in Fig. 1111 
I note that close to the lowest resonance frequency of the 
dispersion at f lja/c ~ 0.6174, there are a large number 
of RSs approaching the pole from smaller frequencies. 
These RSs arise because the refractive index is diverging 
to positive infinity on the low frequency side, allowing for 
a countable infinite number of WG and AWG modes to 
form. <jj is negative for all the resonances in Eq. (1531) . 


Appendix C: Calculation of L 


wavelength X (pm) 
10 1 



RS frequency Re(aco/c) 

FIG. 11: (a) RP as a function of frequency given by Eq. (1531) . 
(b) Resonant frequencies ui n of the homogeneous dielectric 
waveguide with pa — 5 and this RP, forming the basis states 
for the Sellmeier RSE discussed in Sec. IVI El (c) Zoom of (b) 
showing the series of poles on the low frequency side of the sin¬ 
gularity in RP. I show in (b) and (c) states with both incom¬ 
ing and outgoing boundary conditions, i.e. having positive or 
negative imaginary parts of complex resonant frequency, only 
states with outgoing boundary conditions go into the basis. 


Here I derive the form of the L operator occurring 
in the reduced Maxwell’s wave equation for waveguide 
Eq. (JT]). 

For planar systems L is given in Ref.@ as 

= < C1 > 

For effectively two-dimensional systems with one direc¬ 
tion of translational invariance L can be calculated as 
follows. Due to translational invariance in the direction 
z, the direction of propagation I can write the full field 
as, 

E(r) = e ipz ~ iu]t F,(r) = 0 p E(f), (C2) 

where r give the coordinates in the plane perpendicular 
to 2 . 

Therefore using standard vector identities I can write 

VxVxE = VxVx (0 p E) 

= V x (V0 P x E + 0 P V x E) 

= V(9 P (V • E) - E(V • (V0 P )) 

+(E • V)V0 p - (V0 P • V)E 
+V0 P x V x E 
+0 P V x V x E 


I make use of the following identities, 

V0p = ipd p e z (C4) 

V • V0 P = -p 2 0 p (C5) 

V0 p (V • E) = ipO p (S7 ■ E)e z (C6) 

(E • V)V0 p = (e z -E )p 2 d p e z (C7) 

<9 Z E = 0 (C8) 

where e z is the unit vector in the x direction so as to 

simplify Eq. (IC3I) to 

VxVxE = V x V x (E0 P ) (C9) 

= ip9 p e z (y ■ E) + Ep 2 0 p 


+(e z • E )p 2 0 p e z - ip9 p d z E 
+ip9 p e z x V x E + 0 P V x V x E 
= 9 p L{ E) 


(C3) 
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Hence I see that, 


has a Laurent series at q n . Substituting the expression 
Eq. (1D2I) into Eq. (ID3I) gives 


L(E) = ipe z (V ■ E)+ Ep 2 (CIO) 

+(e z • E)p 2 e z + ipe z x V x E 

+V x V x E 


lim (k - q n ) Y] , r ^ = Q„(r, r' 

fc->9n z — K — q m 


so that 


(D4) 


So L is a linear operator independent of the RP and 2 
as required. Therefore the eigenvalue problems derived 
in Sec. IIIII and Sec. m are valid for planar and cylindrical 
wave guides. Also the derivation of L demonstrates I can 
use separation of variables in Eq. (O to arrive at Eq. (IC2I) . 


lim {k - q n ) ^2 ^ = 0 

k->q n k — q m 

m^n 


(D5) 


Substituting the expression Eq. (ID2I) into Eq. m and 
convoluting with an arbitrary finite vector function D(r) 
over a finite volume V we obtain 


Appendix D: Spectral representation of the GF of 
an open system 


E 


—LF n (r) + e(r, k)k 2 F n (r) 
k- q n 


D(r), 


(D 6 ) 


Here I almost exactly repeat my derivations which I 
contributed to Ref. Q using exactly the same method but 
with increased mathematical rigour, in order to prove in 
this section the spectral representation of the Green’s 
function (GF) of a general wave equations. 

The GF of an open waveguide equation system is a 
tensor function Gk which satisfies the outgoing wave BCs 
and the waveguide wave equation Eq. (j4|) with a delta 
function source term, 


where F n (r) = f y Q„(r, r')D(r')dr'. Multiplying by (fc— 
q n ) and taking the limit k —>• q n yields 

lim (t _ gn) y -LF'»(»+mW*W) (D7) 

^ k — q m 

= lim (k — g n )D(r) = 0 . (D 8 ) 

fc— 

From Eq. (ID5I) we can see, 


- LG fc (r, r') + e(r, k)k 2 G fe (r, r') = i<5(r - r') , (Dl) 

In this Appendix the effective permittivity e(r, k ) = 
e w (r)(l +p 2 /k 2 ) with fixed p. Assuming a simple-pole 
structure of the GF inside the scatterer with poles at 
k = q n and taking into account its large-fc vanishing 
asymptotics, the Mittag-Leffler theorem allows us to ex¬ 
press the GF only inside the scatterer as the convergent 


, v- -TF m (f)+ e(r,fc)fc 2 F m (f) 

£%!* ~ fc) A--= 0 • 

777-7=71 

(D9) 

so we can drop terms n 7 ^ m from the summation in 
Eq. (ID 8 I) to give 


lim {k- q n ) 
k—*qn 


-LF n (r) + e(r, fc)fc 2 F n (r) 
k- q n 


= 0 . 


(DIO) 


G fe (r, r') = 

n 


Qn(r,rQ 
k - q n 


or 

(D2) 


- TF„(r) + e(r, q n )q 2 n F„(r) = 0 . 


(Dll) 


It will probably be that Eq. (ID2I) will need to be exper¬ 
imentally verified by comparing scattering predicted by 
the RSE Born approximation with experimental results 
of scattering from well defined scatterers. It may be that 
Eq. (lD2l) is a fundamental law of Physics. My justifi¬ 
cation for this form of the GF is the superposition of 
Lorentzians which make up the scattering profile of res¬ 
onators and the numerical verification of this form of GF 
made in Ref. 

Assuming no degeneracy with the mode n, the defini¬ 
tion of the residue tensor Q„(r, r') at a simple pole of 
the function G*(r, r') is, 

lim [k - g„)G fc (r,r') = Q n (r,r') (D3) 

I have again assumed Gfc(r, r') to be holomorphic in this 
neighbourhood of q n except for at the poles q n so that it 


Due to the convolution with the GF, F„(r) satisfies 
the same outgoing wave BCs. Then, according to Eq. (U), 
F„(r) oc E„(r) and q n = k n , i.e. 

LE„(r) = e(r, k n )k 2 n E„(f), (D12) 

Note that the convolution of the kernel Q n (r, r') with 
different vector functions D(r) can be proportional to one 
and the same vector function E„(r) only if the kernel has 
the form of a product: 

Q n (r, t') = E„(r) ® E„(r')/2A: n , (D13) 

where ® is the dyadic product operator. 

The symmetry in Eq. (ID13I) follows from the reci¬ 
procity theorem, described mathematically by the rela¬ 
tion 

SiG fc (ri,f 2 )s 2 = s 2 G fc (r 2 ,ri)si, (D14) 
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which holds for any two point sources s i t 2 at points r 
emitting at the same frequency. Hence G k (j, r') is sym¬ 
metric. 

In the case of a GF made up of degenerate modes the 
proof of Eq. (ID ldp is modified by making use of orthogo¬ 
nality of the degenerate modes to choose D(r) such that, 


Combining Eq. (IE1D and Eq. (IE6I) yields 


Gfc( r,r') = 

n 


E„(r) 0E„(r') 
2 k(k — k n ) 


(E7) 


Combining Eq. (IE3I) and Eq. (IE6I) leads to the closure re¬ 
lation 


E m (r) ■ D(r) dr = 0 , 


(D15) 


for m ^ n and where state m is degenerate with n. 
Hence I obtain 


G 


k( r,r 


') = 




E»(r) ® E ra (r') 
kn) 


(D16) 


Appendix E: Derivation of sum rule and 
completeness 


—7 p- En(r) ® E„(r') = id(r — r') , (E8) 

n 

which expresses the completeness of the RSs, so that any 
function can be written as a superposition of RSs. If in 
the perturbed system some of the series are not conver¬ 
gent or are instead conditionally convergent then we will 
not arrive at the sum rule and completeness, in which 
case I expect that the RSE Born approximation will still 
give convergence to the exact solution but only if a valid 
spectral Green’s function is used, such as Eq. (IE1I) . 


Here I exactly repeat my derivation of the sum rule 
and completeness of the GF which I made for Ref. |3j but 
in more detail. 

In order to simplify the RSE Born approximation we 
require an appropriate spectral form of the GF which is 
different from the one already proven in Appendix D. To 
obtain this correct form I start with the GF valid inside 
the scatterer only, 


Gfc(r,r') = E 

n 


E»(r) ® En(r') 
2 kn{k kn) 


(El) 


Substituting Eq. (IE1D in 

-LG k (r, r) + (k 2 +p 2 )£ bJ (r)Gk(r, r) = lS(r-r '), (E2) 


gives for k —> oo, 

(k + k n ) E„(r) ® E„(r') 


i(f) £ 


2 k„ 


= i<5(r - r'). (E3) 


since throughout the derivation in this appendix we are 
considering the limit where k —> oc at which e k { r) = 
e(r), i.e. the system is non-dispersive at high frequencies. 

Convoluting Eq. (IE3|) with arbitrary finite functions D 
which are only non-zero inside the resonator, gives 


Kf-lE 


(■k + k n ) E n (r) ® E„(r') 
2 k n 


D(r , )dr / = D(r). 

(E4) 

and assuming the series are convergent as k —> oo we get, 
E„(r) <g> E„(r') 




2 kn 


D(r')dr' = 0. (E5) 


and since this is true for all D it follows that 


E„(r) ® E„(r') _ ^ 


2 kn 


(E6) 


Appendix F: Derivation of the correct normalisation 
for effectively 2D waveguides 


Here I almost exactly repeat my derivations which I 
contributed to Ref.@] using exactly the same method as 
for Maxwell’s equations in order to prove in this sec¬ 
tion that the spectral representation, only valid inside 
the waveguide (similar to Ref.[3] for nano-particles), rig¬ 
orously derived in Appendix C, D, and E 


G fe (r,r') = 

n 


E„(r) ® E ra (r') 
2 k(k — k n ) 


(FI) 


leads to the RS normalization condition Eq. (IF8I) for gen¬ 
eral waveguide equations. To do so, I consider an an¬ 
alytic continuation E(r, k) of the vector wave function 
E„(r) around the point k = k n in the complex fc-plane 
( k n is the wavenumber of the given RS). I select the ana¬ 
lytic continuation such that it satisfies the outgoing wave 
boundary condition and my waveguide equation (taking 
c= 1) 

—AE(r, fc)+e(r,w)(/c 2 +p 2 )E(r, k) = (k 2 —k 2 )cr( r) (F2) 


with an arbitrary source term. 

The source <x(r) has to be zero outside the cross-section 
of the inhomogeneity of £ w (r) for the field E(r, k) to sat¬ 
isfy the outgoing wave boundary condition. It also has to 
be non-zero somewhere inside that cross-section, as oth¬ 
erwise E(r, k) would be identical to E n (r). It is further 
require that <r(r) is normalized according to 


E„(r) • er(r) dr = 1 + 5 kn ,o ■ 


(F3) 


The integral in Eq. (IF3I) is taken over an arbitrary cross- 
section A which includes all system inhomogeneity of 
e w (r). I do not derive Eq. (IF3I) . it is simply a conve¬ 
nient condition to put on the otherwise arbitrary spatial 
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dependence of the source which lies only inside the waveg¬ 
uide scatterer. If I had made the condition anything else 
(which in my earliest proof I did) then algebra and can¬ 
cellation would lead us back to the same result, however 
with more mathematical complexity and operations. The 
condition make the mathematics easier because it causes 
E —> E„ exactly and without proportionality constant 
appearing, equation (IF3I) ensures that the analytic con¬ 
tinuation reproduces E„(r) in the limit k —> k n . Solving 
Eq. (1F2I) with the help of the GF and using its spectral 
representation Eq. (IF1I) . I find: 

E(r, k) = f G fc (r,r')(/c 2 - k 2 n )a(v')dv' 

J A 

(F4) 

and then, using Eq. (IF3I) . obtain 

lim E(r, k) = E„(r), 

k — yk n 

for any r inside the system. Outside the system, the 
analytic continuation E(r, k) is defined as a solution of 
the waveguide equation wave equation in free space. This 
solution is connected to the field inside the system [given 
by Eq. (IF4j) ] through the boundary conditions. Note that 
in the case of degenerate modes, k m = k n for m ^ n, 
the current er(r) has to be chosen in such a way that it 
satisfies Eq. (IF3I) and, additionally, 


[ E m (r) • cr(r) dr = 0, 
J A 


in order that the degenerate modes can be normalised 
separately. 

I now consider the integral 


and then finally with La the circumference of A 1 


l + <V-„,o — / E„(r) 


duj 2 e(r,uj) 


du 2 


E„(f)dr(F 8 ) 


, E„ • VE — E ■ VE„ 
lim ® -—--- aL . 


k—tk n 


La 


U 2 _ jU2 


I have also assumed a simpler form for the waveguide 
equation in free space, see Appendix G so that I can 
use the vector and scalar Green’s theorem to simplify 
En.(lF7lt. Various schemes exist to evaluate the line 
inte gral limit in Eq. (IF8I) such as analytic methods in 
Ref. 11] or numerically extending the closed arc into a 
non-reflecting, absorbing, perfectly matched layer where 
it vanishes. 

Due to the over completeness of the basis demonstrated 
in Appendix E we do not have the relation we find for 
Hermitian systems 




i(*0 — 


£u 


>(*) 


' allspace 


E m (r) • E„(r) dr = 0 , 
(F9) 


for all to ^ n. The basis states of non-Hermitian sys¬ 
tems are not all mutually orthogonal although some are 
orthogonal to one another. Hence we cannot derive for 
non-Hermitian systems the normalization relation, which 
is well known for Hermitian systems and follows from 
Eq. (|F9j) . 


f d(u; 2 s^(f)) 

' allspace d{uj ) 


E„(r) • E„(r) dr = 1 + S kn ,o , 

l n 

(F10) 


and, hence, Eq. (IF8I) and Eq. (IF10D are different. 


Appendix G: RSE Born Approximation for 
effectively 2D waveguides 


In(k) = 


fA (E • TE n — E„ ■ TE)dr 

7^2 _ U2 

h, 


(F5) 


In effectively two dimensions Eq. ([2]) becomes in free 
space (taking c = 1), 


and evaluate it by using the waveguide equations Eqs. Q 
and (IF2I) for E n and E, respectively, and the source term 
normalization Eq. (IF3I) : 


T /a(E ■ w 2 £a, n E„ — E„ • w 2 e w E)dr 

4(fc) = —- Y2 - U2 -H + 4.,0, 

K K n 

f F6 ) 

where I assume that is a a (real) symmetric matrix or 
a scalar so that (defined to be in this case) E ■ e Wn E„ = 
E„ • e Uii E, and so I obtain also by commutation of E and 
E„ and simple calculus in the integral shown Eq. (IF6I) 
the dispersion factor and normalization becomes, 


1 + dfc o — 


En(f) • 


d(u} 2 £ u (r)) 


lim 


d(uj 2 ) 

E n (f)LE(f,k)-E(r,k)LE n (r) 


k—tk- 


n J A 


VI _ U2 

ib Ibm 


E n (f)df (F7) 

dr 


V 2 - 


1_ 

c 2 dt 2 


d P E(f)=0, 


(Gl) 


However V 2 d p = —p 2 9 p , therefore 

[V 2 + fc 2 ] E = 0, (G2) 

Hence the free space GF equation is, 

[V 2 + fc 2 ] G-{ s (r, y') = id(r — r') , (G3) 

which has the solution, also assuming p >> p ', 


G{ s (f,f') = -lfF 0 (fc|r-f'|)i « !1 (G4) 


The systems associated with G;[ s and G k of Eq. m are re¬ 
lated by the Dyson Equations perturbingback and forth 
with Ae^(r) = e u (r) — i similar to Ref. ( 4 ], 
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G fc (r,r") = G{ s (r,r") 


k v*!* ) (G5) 

- W 2 / G{ s (r,f"')A£ I ,(r'")G fc (r'",r")dr'", 


G fc (r"',r") = G[\v'\v") (G6) 


-c, 2 / G k (i"',i')Ae u (i')G{'(T',i'')di f 


Combining Eq. (IG5I) and Eq. (IG6I) it is obtained similar 
to Ref-HI 

G fe (r,f") = G{ s (r,f") 

- w 2 j G{ s (f, r , )A£ w (r')G{ s (r , I r")dr' 

+ w 4 J j G{ s (r, r , )Ae w (r / )G fc (r', r'") 

xAe u (r"')G{'(r"'r ")dr"'dr'. (G7) 

In order to improve the numerical performance further 
I make a final few steps as in the original Born approx¬ 
imation [I], I define unit vector r such that r = pr and 
k, = hr. Then for p » p\ 


k\v - r'| = kp - k s • v' + O ( — 


P- 


(G8) 


Therefore substituting Eq. © and Eq. (1G4I) in to 
Eq. (IG7I) and using Eq. (IG8I) because both r, r" are far 
from the scatterer I arrive at the RSE Born approxima¬ 
tion 


G k (r,r") = -^H 0 (k\r-r"\)l 


— UJ 


j Ae u (r')dr' 


i Q 2 e ifc(p+p" 

rVpp* 

A Q2 e ik( P+P ") ^ A„(k s )® A„(-k s ") ^ 

+ W TTV^ 2 kik-kn) • (G9) 

or using Eq. (EH) instead 


I note that the fast Fourier transform method is avail¬ 
able. The first two terms in Eq. (1G9I) correspond to the 
standard Born approximation, the final summation term 
corresponds to the RSE correction to the Born approxi¬ 
mation. 

A simple corollary of this theory is as follows, I can see 
from the arguments just stated that from Eq.(G5) if r" 
is inside the resonator and p » p" then 

G fc (i,i'0 = —>o(fc|r-r"|)i 


— UJ 


sfp^ 

or using Eq. EH instead 


A n (k s ) g) E„(r") 

A- 2 k(k-kn) ’ ’ 


G fc (f,f") = - 7 ^o(fc|r-r"|)i 


— UJ 


J2 An ^ s) n® En{i ; n ' . (Gi3) 


y/pn ^ 2 kn{k - k n ) 


similarly from Eq.(G6) if r is inside the resonator and 
p" » p then 

G k (r,r") = - l -H 0 (k\r-r"\)i 

2 Qe ik »" ^ E n (r)®A n (-k a ") ^ 

2k{k _ kn) - ( G14 ) 

or using Eq. (1E1I) instead 

G fe (r,f") = -^o(fc|r-f"|)i 

2 Qe ikp " r E »(f)®A„(-k;j 

_ w - KT77T, - — -> (GI5J 




n 


2 k n {k k n ) 


other permutations are possible. 

It may be possible to use these formulas to develope 
a hollow glass cable with light and fluid passing along 
it, the escaping light being measured in the far field for 
analysis with the RSE Born approximation for the deter¬ 
mination of the content of the fluid via inverse emission 
problem. To elaborate on this point, it was shown in 
Ref.[|| that 


Gfc(r, i") = ~^H 0 (k\r — f"|)l 


i Q2 e ifc(p+p") 


+ A 


t Vpp 77 


y~~ A ra (k^) ® A„(—k/') (G10) 


iQ 2 e ik(p+p") A „(k s ) g A „(—k s ") 


2 kn{k kn) 


The vector A n is defined as a Fourier transform of the 
RSs, 


A„(k s ) — 


, / k 


Ae w (r')E n (r')dr' 


(Gil) 


lim E„(r) oc A„(ffc n )^— j- . (G16) 

p—YOO p / 

However close to a sharp resonance, when [k ~ 
p is fixed], scattering is dominated by a single resonance 
and so the scattered E-field E scattered is approximately 

* p ik P 

lim E scattered (r, k) m CA n (rk) — . (G17) 

p—too p*-/* 

Hence A„ (f k) can be partial inverse Fourier transformed 
with respect to angle to find information about the in¬ 
ternal structure of the cable, C is some constant. The 
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same approach can be taken in three dimensions to aid 
the determination of resonator structure. 

The reason Eq. (IGlfijl takes the form it does can be seen 
more clearly by substituting the GF of Eq. (IG12D into the 
first line of Eq. (IF 41) and taking the limit k —» k n , while 
following the arguments in that section that, 

lim E(r, k) = E„(r). 

k y k n 

The 3D equivalent i! of Eq. (IG16I) one day might 
be valuable information for solving the inverse emis¬ 
sion problem from such things as black hole gravita¬ 
tional wave emitters, i.e. calculations of the poten¬ 


tial from the emission (decay) via fast inverse Fourier 
transform methods upon the set of A n (r k n ), particu¬ 
larly if the potentials of interest are rotating about a 
fixed axis so we know their orientation to some extent 
such as occurs for decaying magnetic nuclei as part of a 
non-magnetic crystaline compound placed inside a NMR 
(nuclear-magnetic-resonance) machine. Because k n are 
discrete values A„(r k n ) only give angular information 
when inverse Fourier transformed and so the inverse 
Fourier methods might have to be used self-consistently 
in conjuncture with the RSE perturbation theory and the 
values of k n . This is a highly speculative aside and might 
be a possible topic for future research. 
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